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CHAPTER I 


INTRODUCTION 

A. BACKGROUND 

Satellite' communication is suitable for a wide, 1 variety of 
applications. Its prime use is to provide long-distance communications, 
and this application is continuing to grow fast. 

One of the reasons why satellite communication has become so 
popular is the unique synchronous property of the geostationary 
satellite orbit (GSO). The GSO can be described as an imaginary circle 
in the equatorial plane, approximately 22,300 miles (36,000 km) above 
the surface of the Earth. A satellite placed there orbits the Earth at 
the same angular velocity as the Earth rotates on its axis, so that it 
appears to be stationary to an observer on the Earth. Although there 
are some important satellite applications using non-geostationary 
orbits, without doubt the existence of the GSO has made satellite 
communication so significant. 

Satellites cannot be placed arbitrarily close to one another along 
the GSO, not primarily because of danger of possible collision, but 
because of the resulting electromagnetic interference. Allowable 
interference levels dictate the minimum possible spacings between two 
satellites using the same frequency band. Therefore, the GSO can 
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accommodate only a finite number of satellites using any given frequency 
band. It has been recognized that the GSO is a limited natural resource 
[1], As use of the GSO becomes more intensi ve, the effective 
utilization of this natural resource has become an important issue. 

Technical analyses on satellite communication services are 
performed by the International Radio Consultative Committee (CCIR) of 
the Internationl Telecommunication Union (ITU), using the work of CCIR 
study groups as input. This study utilizes the CCIR recommendations and 
reports as sources for technical parameters, such as antenna reference 
patterns. 

The purpose of this study is to develop computational methods and 
to perform engineering calculations which will provide a basic 
understanding of satellite orbital locations. More specifically, 
minimum satellite separations are calculated which satisfy a given 
carrier-to-interference power ratio (C/I) for the Fixed-Satellite 
Service (FSS) on a single-entry basis. The results are presented in 
sets of contour curves as well as computer codes. 

B. SOME PRELIMINARIES 

A satellite communications link consists of an up link from an 
earth station to a satellite and a down link from the satellite to a 
destination earth station. Interference can occur on each of these half 
links. On the up link, interference may be due to the signals 
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from other up links, from other', satellites,' pr 'from^s^tioris, in .the 
terrestrial service. In this study, it wi 1 V be assumed that the up -link 
signals under consideration will be in di ffereht -frequency bands from 
all down-link, inter-satellite, and terrestrial -service signals, or that 
separate coordination procedures will be used, sp that in up-link 
calculations only the interference between up links .needs to be 


considered. 



under consideration will be in frequency bands di ffereht' from those 
assigned to up links, inter-satellites, and terrestrial services, or 
that these services will be coordinated, so that in. down-1 ink 
calculations only the interference between down links need be 

„ , 'X 

considered. Therefore, each half link will be treated separately. The 
situation is shown schematically in Figure 1.1 with the geometric 

• ■ •• *■ . - *s • 

parameters of interest. •' \ 

As shown in Figure 1.1, only two satellite networks will be dealt 
with in this study, i.e., the analysis is on a single sentry basis. One 
is viewed as a wanted network, and the other as an. -interfering network. 
Earth station antennas are assumed to be pointed directly at their 
respective satellite antennas as required by good system design. Also, 
pointing and station keeping errors are ignored for simplicity, and only 
circular-beam antennas are considered here. 

In actual system analysis, many interfering satellites or 


interfering earth stations are involved (aggregate. interference) rather 


than just a single interfering source (single-entry\interference) . In 
order to obtain an estimate of multi -interference source situations. 



WANTED SATELLITE (SW ) INTERFERING SATELLITE (SI) 



EARTH STATION EARTH STATION ' EARTH STATION EARTH STATION 

FOR WANTED FOR WANTED ' FOR INTERFERING FOR INTERFERING- 

SATELLITE SATELLITE SATELLITE SATELLITE 

( EWR ) ( EWT ) (EIR) ( EIT) 

WANTED NETWORK INTERFERING NETWORK 

WANTED SIGNAL PATHS 

INTERFERING SIGNAL PATHS 


Figure 1.1 Interference geometry between satellite networks. 
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the carrier-to-interference power ratio (C/I), which will be used here 
as an interference criterion, should be raised by a suitable factor, 
such as 5 dB. For example, if the required C/I ratio on an aggregate 
basis is 30 dB, a single-entry C/I ratio of 35 dB due to one of the 
worst interference sources will often be satisfactory [2], 



CHAPTER II 

DOWN-LINK INTERFERENCE 

A. INTERFERENCE GEOMETRY 

The single-entry C/I ratio will be derived first for down links. 

The interference geometry between the down-link parts of two FSS 
networks is shown in Figures 2.1 and 2.2, where ^ 3 , and i |>5 denote 

angles between the main-beam axis and the direction considered (off-axis 
angles), and lowercase letters a, c, e, and g denote distances. The 
wanted signal is received from satellite SW, while the interfering 
signal arrives from satellite SI. The calculation of the C/I ratio is 
performed at earth station EWR, and the resulting equation is used to 
determine the allowable satellite spacing. The ground aim point of 
satellite Si is designated by ElR; an earth station of the interfering 
network may (or may not) be located there. Satellite SW is not 
necessarily pointed directly at earth station EWR; however, earth 
station EWR is located within a service area of satellite SW. 

B. CARRIER AND INTERFERENCE POWERS 

It is first necessary to determine the carrier signal power and the 
interference signal power at earth station EWR. Signal power levels can 
be obtained by application of the Friis transmission formula [3] 



WANTED SIGNAL PATHS 

--INTERFERING SIGNAL PATHS 


Figure 2.1 Interference geometry between down-link networks 




GEOSTATIONARY 
SATELLITE ORBIT 



Figure 2.2 Alternative presentation of interference geometry between down-link 
networks (see Figure 2 . 1 ). Earth radius exaggerated for clarity. 









( 2 . 1 ) 


PR = 


P T G I G R: X > Y , 
(4 tt r)^ 


where P denotes power, G gain, X wavelength, and r distance. The 
subscripts T and R indicate transmit and receive, respectively. 

The carrier signal power and the interference signal power 
l£WR | ~ecei ved by earth station EWR are given by 


'EWR 


P SWT {G SWT D SWT^5 ,G SWT^ }G EWR X W 

( 4 ) 2 


and 


( 2 . 2 ) 


EWR 


P SIT {G S1T D SIT^2’ G SIT^ }{G 


(4ttc)‘ 


EWR°EWR 


(to,G 


EWR^ }X I 


(2.3) 


where the multiple subscript notation adheres to the following 
convention: the first letter identifies the location, S for satellite 

and E for earth station, of the parameters in question; the second 
letter indicates wanted (W) or interfering (I) network; and the third 
letter indicates transmit (T) or receive (R). Thus, the terms in 
Equations (2.2) and (2.3) are interpreted as follows: 

P SWT = power transmitted from the wanted satellite 

p S IT = power transmitted from the interfering satellite 

g SWT = transmitting gain of the wanted satellite antenna 
in the direction of the main beam axis 

Gs IT = transmitting gain of the interfering satellite 
antenna in the direction of the main beam axis 

G EWR = receiving gain of the wanted satellite antenna 
in the direction of the main beam axis 
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= angle between the axis of the interfering 
satellite antenna and the direction of the-; 
wanted earth station , - - ■ 

1J/3 = angle between the axis of the wanted earth 
station antenna and the direction of the . 
interfering satellite 

ij>5 = angle between the axis of the wanted satellite 
antenna and the direction of the wanted earth 
station 

d SIT( ^GsiT) = relative gain of the interfering 
satellite antenna in the direction 
of the wanted earth station 

d EWr( 'KJ’Gewr) = relative gain of the wanted earth 

station antenna in the direction of 
the interfering satellite 

d SWT( ^»^SWT) = relative gain of the wanted satellite 
antenna in the direction of the wanted 
earth station 

= wavelength of the wanted signal 

Xl = wavelength of the interfering signal 

a = distance between the wanted satellite and the wanted 
earth station 

c = distance between the interfering satellite and the 

wanted earth station 


The units Watts, degrees, and meters 

C. REFERENCE ANTENNA PATTERNS 

The relative gain DsiT( 1>2 »GsIT) 
angle <j£ and the gain G$u indicated 
from the satellite antenna reference 
the relative gain Dewr(i1>3,Gewr) is a 


are implied above. 


is a function of the off-axis 
in the argument and is determined 
pattern shown in Figure 2.3, while 
function of ^3, and G^wR ar, d "> s 
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determined from the earth station antenna reference pattern in Figure 
2.4. These reference patterns are adopted from CCIR Report 558-2 [4] 
and from CCIR Report 391-4 [5], respectively. Note that the reference 
patterns in Figures 2.3 and 2.4 are expressed in decibels. 

Consequently, Dsu, G$it, DgwR. and G^^R in these figures are all in dB, 
not ratios. The notation dB is placed before the equation number to 
indicte a decibel form throughout this article. 

The half -power beamwidth i|q in Figure 2.3 can be calculated from 
the gain Gsjj by the relationship 


% 



27000 
G S I T 


degrees 


(2.4) 


as suggested in CCIR Report 215-4 [6]. 

The antenna di ameter-to-wavelength ratio d/X in Figure 2.4 is 
related to gain G^r as 


^EWR ~ eiT » (2.5) 

where the aperture efficiency e is assumed to be 55% for circular 
parabolic reflector antennas [7]. 
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RELATIVE ANTENNA GAIN 



°SIT = 

-12 ( \p2/ 4o) 2 dB 

for 

0° < 1 P 2 < 1 . 3 \}>q 

(S-l) 

°SIT ■ 

-20 dB 

- - . — 

for 

1. 3\po < 4*2 < 3.154b 

(S-2) 

d SIT = 

-7.5 - 

25 1o 9 lfl 

( / 4»o) dB 






for 

3 . 1 5 \pg < 4>2 < 4*1 

(S-3) 

°SIT = 

-Gsit ■ 

- 10 dB 

for 

4i < 4£ 

(S-4) 


where 4 q : half -power beamwidth 

4>j : value of ^ w ^ en &SIT region (S-3) is equal to 

-Gsit - 10 dB 


Figure 2.3 FSS satellite antenna reference pattern. This reference 
pattern is taken from CCIR Report 558-2 [4]. 
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RELATIVE ANTENNA GAIN 



0 \j/ r 48 ° 


OFF-AXIS ANGLE, \f/ (DEGREES) 


'ewr 

= -2.5 x 10 " 3 (d/x ) 2 *, 2 dB 

for 

0 ° 

< ^3 < 


(E-l) 

EWR 

= D x dB 

for 

4> 

< ^3 < 

V 

( E — 2 ) 

EWR 

= 32 - 25 log iQ ( i' 3 ) - ®EWR 

for 

ib 

r r 

< *3 < 

48° 

( E — 3 ) 

'ewr 

= -G ewr - 10 dB 

for 

O 

00 

1 < ^3 


( E— 4 ) 


where 

d : antenna diameter. 

X : wavelength 

Di : relative gain of the first sidelobe = 2 +15 log 10 (d/x) - Gewr 
= 20 /-D i / (d/x) degrees 
’V = 15.85 (d/X)"0*6 degrees 

Figure 2.4 FSS earth station antenna reference pattern. This reference 
pattern is taken from CCIR Report 391-4 [5]. 
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D. CARRIER-TO-INTERFERENCE POWER RATIO 


The single-entry C/I ratio at earth station EWR is obtained by 
dividing Equation (2.2) by Equation (2.3) as 

P SWT G SWT D SWT^5 ,G SWT^ c 2 

( / ^ EWR P SIT G SIT D SIT^2 ,6 SIT^ D EWR^3 ,G EWR^ a 2 ' 
where it is assumed that the carrier signal and the interference signal 
use the same frequency, i.e., Ay = Xj. 

A small notational simplification is available in Equation (2.6) by 
denoting 

E SWT = P SWT G SWT ( 2 * 7 ) 

and 

E SIT = P SIT G SIT » ( 2 * 8 ) 


which are the effective isotropically radiated powers (EIRP) of the 
satellite transmitters in the direction of their respective beam 
maxima. Equation' (2.6) now becomes' 


(C/I) 


’SWT 


D SWT^5, G SWt) 


EWR 


E SIT D SIT^2’ G SIT^ D EWR^3’ G EWR^ a 2 


(2.9) 


For reasonable antenna patterns, the minimum satellite spacing 
angles calculated by use of Equation (2.9) turn out to be relatively 
small, and consequently the distances a and c are nearly equal, i.e., 
c2/a2 is in the range of 1 ± 0.028 for 5° or smaller satellite 
separation. Note that Figures 1.1, 2.1, and 2.2 are not drawn to scale. 
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The actual ratio between the radius of the Earth and the radius of the 
geostationary satellite orbit is 1 : 6.6. Thus, by assuming c 2 /a 2 = 1, 
Equation (2.9) can be simplified to 
E, 


(C/I) 


'SWT 


D SVIT^*5 ,G SWT^ 


EWR E $it d SI t^2 ,6 SIT^ D EWR^3 ,G EWR^ 


( 2 . 10 ) 


E. UNIVERSAL SYSTEM PARAMETER 


Since the objective here is to find the minimum possible spacings 
between satellites and to present the results in a form which is helpful 
for system design, it is useful to reduce the number of parameters as 
much as possible in Equation (2.10) by combining some of them into a new 
system parameter. Defining such a parameter Rdn as 

E SIT 1 (2.11) 

DN ( C/I )EWR E swt D swt (%»S swt ) 

simplifies Equation (2.10) to 

= l _ ( 2 . 12 ) 

DN d sit^2 ,G SIT^ D EWR^3 ,G EWR^ 

Equations (2.11) and (2.12) can also be expressed in a decibel form as; 

r dn = ( c/i )ewr ‘ E SWT + E SIT “ D SWT^5 ,G SWT^ dB ( 2 * 13 ) 

and 

R DN = " D SIT^2’ G SIT^ " D EWR^3 ,G EWR^ dB ’ (2.14) 

respectively. The decibel form equations will be used in the 
subsequent discussion. 
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For each network system, the required (C/ 1 )ewr> EsWT> ancl ^SIT 1 n 
Equation (2.13) have given values, while DsWT depends on ^5 and G$wj. 
However, Oswt varies only from 0 to -3 dB, i.e., D$WT is ignored (0 dB) 
for EWR located at the SW aim point and is assumed to be -3 dB for EWR 
at the edge of the SW coverage area. Thus, D$wt in Equation (2.13) can 
be estimated between 0 and -3 dB, depending on how closely EWR is 
located to the SW aim point. 

Equation (2.14) can apply not only to the co-channel interference 
case, but also to the adjacent -channel interference case by obtaining 
the required C/I ratio for the adjacent-channel interference. An 
example is provided in Figure 2.5, which gives the adjacent channel 
protection ratio relative to the co-channel protection ratio as a 
function of the frequency offset. For example, suppose two channels 
are separated by 20 MHz, i.e., the center frequencies of the channels 
differ by 20 MHz, then the required C/I ratio is 8 dB below that for the 
co-channel interference case. Thus, Rqn calculated from Equation (2.13) 
is also 8 dB bel ow that for the co-channel case. 

The parameter Rqn will be used as a criterion reference level and 
called a universal system parameter in the subsequent discussion. It 
can also be interpreted physically as the total antenna discrimination 
from Equation (2.14). In application, the left side of Equation (2.14) 
is easily evaluated from Equation (2.13), while the right side is 
evaluated by computer and is presented in the form of a graph, i.e., the 
contour curve of ^3 against ^2 a t each value of Rqn can plotted for 
fixed values of G$^t and Ge^/R by using Equation (2.14). The computer 
code for this implementation is given in Appendix A. 
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Figure 2.5 Protection ratio relative to co-channel values [8], 
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F. GRAPHICAL PRESENTATION 


The result from Equation (2.14) is presented graphically in Figure 
2 . 6 , where the minimum allowable satellite spacing angle 1^3 is shown as 
a function of the off-axis angle ife as seen from satellite SI for Ron = 
30 dB, 35 dB, and 40 dB. The system parameter values G$it = 50 dB and 
%WR = 50 dB, corresponding to a half-power beamwidth of approximately 
0.5°, are used for the calculation in Figure 2.6. 

The off-axis angle ^ is the angle between the main-beam axis of 
the SI antenna (toward EIR) and the direction of EWR, and can be 
interpreted as an earth station separation between EWR and EIR as seen 
from SI. It can be calculated from the two earth station locations and 
the orbital satellite location in terms of longitude and latitude by 
the Law of Cosines, i.e., 



cos 



+ e 2 - 


2 -ce 



> 


(2.15) 


where the distances c, e, and g, defined as in Figures 2.1 and 2.2, are 
obtained in the following way. Referring to Figure 2.7, the locations 
of EWR and EIR are represented in the rectangular coordinate form as 


(rcos 9 ewRCOS<}>ewr, rcos 0EWR S1 n <I>EWR , rs 1 ' n 9 EWR) 
and 

(rcos 9 Eir cos<})eir, rcos QElRSin^iR, rsin9£iR), 


respectively, and the location of SI is represented as 



MIN. SATELLITE SPACING, 



Figure 2.6 Minimum allowable satellite spacings as function of the 
earth station separation seen from the interfering 
satellite, based on down-link considerations. 

G SIT = g EWR = 50 dB. 
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X EWR = rcos 0EWRCOS 4>EWR 
y EWR = rcosQEWRSin^WR 
Z EWR = rsin0EWR 
x Si = Rcos $si 
y SI = Rsi n <f>si 
Z SI = 0 

Figure 2.7 Relationship between earth station and satellite locations 
in terms of latitude and longitude and their rectangular 
coordi nates . 
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(Rcos 4>si . Rs i n 4>s I » 0) » 


where R indicates the radius of the geostatioary satellite orbit (4.2152 
x 107 m), r indicates the earth radius (6.371 x 10^ m), and 8 and <J> 
denote latitude and longitude, respectively. Note that latitude 8 is 
negative in the southern hemisphere, and that longitude <f> is negative in 
the western hemisphere. The subscript notation follows the same 
convention as before. Thus, the distances c, e, and g are calculated 
from 


R + r - 2Rrcos ( -<j> E y R )cos 8 E ^ R , 

(2.16) 

2 2 

= R + r - 2Rrcos( <j) SI -(|) EIR )cos e EIR , 

(2.17) 

2r {1 “ COS 8^ j R COS ( JR ) 

- sine EWR sine EIR } 

(2.18) 


It is evident that the calculations required to obtain the exact ■[>}';- 
off-axis angle i |>2 are numerous. A simple approximation is available 
for small values of ^ by assuming that c » e. The distance c or e is 
expressed as 5.62r for a satellite operated at a high (90°) elevation 
angle and 6.54r for a satellite operated at a low (0°) elevation angle. 
Thus, 

9/2 

s1n( v2> -iifr-r • 
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g 


radi ans 


(2.19) 


i .e., 

'h “ 5.62r 

for high elevation angle. Similarly, one obtains 
9 

\|>2 ” " 6 ~ 54 r radians ( 2 . 20 ) 

for low elevation angle. Given the distance g, Equations (2.19) and 
(2.20) provide an upper bound of ^2 anc * a lower bound of ^2* 
respectively, and the bounds are quite narrow for small ^ 2 * 

The off-axis angle ^3 is the topocentric (as seen from a point on 
the Earth) satellite spacing angle and not the geocentric (measured from 
the center of the Earth) satellite spacing angle. In practice, we are 
interested in small satellite spacings, and under this condition, these 
two angles are nearly equal. Also, the topocentric angle is always 
greater than the geocentric angle, and hence the minimum satellite 
spacing angle % in- Figure. 2 . 6 is always conservative In thi s sense. — , 
However, the difference between the two angles is not always 
negligible. Especially when satellites are near zenith, and an 
observation point on the Earth is near the equator, the error in 
substituting the topocentric angle for the geocentric angle may 
approach 15%. An approximate but accurate conversion formula from the 
topocentric angle % to the geocentric angle ♦g is given by 

= / 1.023 - 0.302 cos( $M-<t£)cos 0£ radians , (2.21) 

where 4 ^ is the longitude of the midpoint between two satellites, and 0 £ 
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and ‘t’E are the latitude and longitude of an observation point on the 
Earth, respectively. The derivation of Equation (2.21) is shown in 
Appendix B. 

Consider now Figure 2.6. For Ron = 35 dB, the required satellite 
spacing is 4.8° when 4^ = 0°, i.e., when EWR and EIR are collocated. 

The required satellite spacing is reduced to 2.5° when ^ = 0.4°, i.e., 
when EWR and EIR are separated. The condition t 2 = 0° implies that the 
satellite antenna beams are overlapped (overlapping service areas), and 
^2 > 0° implies that the antenna beams are restricted (assigning limited 
service areas). Thus, in the case of ^ = 0°, the EWR antenna is 
illuminated by the main beam of the SI antenna, and discrimination is 
achieved only by the EWR antenna in the near-si delobe region, thus 
resulting in the worst satellite spacing, whereas, in the case of ^2 > 
0°, the beam maximum of the SI antenna does not face the EWR antenna 
directly, and less discrimination by the EWR antenna is required, thus 
allowing closer satellite spacing. Thus, the concept of assigning 
limited service areas is crucial for the effective use of the orbital 
arcs. This concept seems to have been enunciated first by Kiebler [9]. 


G. UNIVERSAL CURVES 

Figure 2.6 is shown for the particular values G$it = GgwR = 50 dB. 
The assignment of specific values to these variables in Equation (2.14) 
allows the plotting of a family of curves with Rqn as a parameter. A 
similar family of curves can be constructed for each set of gain values. 
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but an unwieldy number of curves result. Instead, the curves will be 

reduced to a more universal form. * 

Since the FSS satellite antenna reference pattern of Figure 2.3 is 
represented in terms of ^/^O* where the half -power beamwidth ^ is a 
function of G$jt as shown in Equation (2.4), the abscissa variable in 
Figure 2.6 can be chosen as ^ 2/^0 instead of ^ 2 » thus no longer 
requiring the specification of the value Gg u for plotting the curves. 

On the other hand, the FSS earth station antenna reference pattern 
of Figure 2.4 is represented as a function of (d/x)^ in the region 
labeled (E-l), i.e., 

D EWr(* 3 ’W = -2.5 x 10- 3{ (d/A)V 2 dB * ( 2 . 22 ) 

, 0.4 

and as a function of ^Gewr in the region labeled (E-3), i.e., 

^EWR^ ^3’^EWR^ = 32 - 25 log i o ( ^ 3 ) ~ 10 logjo(^EWR^ 

. : = 32 - 25 log ! o( ^ 3 g ewr) dB » - - - . ( 2 * 23 ) 

in which G^wr is not in dB but is a ratio. Note that G^wr is a function 

of d/Xas shown in Equation (2.5). Thus, the ordinate variable in 

Figure 2.6 can be transformed into (d/X )^3 for the portion of the curves 

,0.4 

contributed by region (E-l), and into ^Gewr for the portion of the 
curves contributed by region (E-3). At first sight, it seems that four 
sets of curves (four different ordinates) are necessary in order to 
present ^3 because there are four regions (E-l) through (E-4) in the 
reference pattern of Figure 2.4. Fortunately, the (E-4) portion is not 



necessary since the ^3 values of interest are the minimum satellite 

spacing angles, which never exceed 48° for any practical situation. 

Furthermore, the curves in Figure 2.6 merely have a discontinuity, as 

shown by the dotted lines in Figure 2.8, when the value falls into 

the region (E-2). Therefore, two sets of universal curves, 

corresponding to the regions (E-l) and (E-3), suffice to present the 

minimum allowable satellite spacing angle ^ 3 . The set of universal 

curves contributed by the region (E-3) is shown in Figure 2.9 with the 
0.4 

ordinate as ^ewr, where the values of ^3 obtained from Figure 2.9 must 
satisfy 

< 1> 3 > t r = 15.85 (d/xf °* 6 

= 26.3 g;2r 3 degrees , (2.24) 


in which G £WR is a ratio. The other set of universal curves, 
contributed by the region (E-l), is shown in Figure 2.10 with the 
ordinate as (d/x)^, where the values of ^3 obtained from Figure 2.10 
must satisfy 


*3 < *m = 20 ^G EWR - 2 - 15 log 1Q (d/X) / (d/X) 


= 20 Y 5 - 35 + 3_log 10 (d/x) / (d/X) degrees.. (2.25) 

The dashed lines in Figures 2.9 and 2.10 indicate where to switch from 
one set of curves to the other for G^ R = 50 dB and 60 dB. The value 
of <|>2 has to be obtained from the area above the dashed lines in Figure 
2,9, or from the area below the dashed lines in Figure 2.10. 
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Figure 2.8 Discontinuity of contour curves in Figure 2.6. 

Both sets of universal curves are based on Equation (2.14). In 
order to construct the first part of the universal curves. Equation 
(2.23) is substituted for D^ WR in Equation (2.14), and (d/x)^ is 
calculated for each rvalue of by means of an iterative search. 

Similarly, in order to construct the second part of the universal 
curves. Equation (2.22) is substituted for D £WR in Equation (2.14), and 
^3 G EWR 1S ca ^ cu ^ ate d for each value of The com P uter codes used 

for producing Figures 2.9 and 2.10 are given in Appendices C and D, 
respectively. 

An example of how to use these curves will be given in the 
following. Given system parameters such as a required C/ I ratio and an 
EIRP ratio, the universal system parameter can be computed from 
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Figure 2.9 First part of universal curves for minimum allowable 
satellite spacing angle ^ as function of normalized 
off-axis angle based on down -link consideration. 

<|>3 found from this figure must satisfy Equation (2.24). 
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Figure 2.10 Second part of universal curves for minimum allowable 
satellite spacing angle ^ as Function of normalized 
off-axis angle based on down-link considerations. 

^3 found from this figure must satisfy Equation (2.25). 
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Equation (2.13). The off-axis angle ^ 1S obtained as explained in the 

previous section F, and the half -power beamwidth ^ is given in each 

system or is obtained from Equation (2.4). Suppose R^ = 35 dB, = 

50 dB, and ^ = 0.15°, then = 0.3. Figure 2.10 cannot be applied 

to obtain ^ in this case because the value of (d/x)^ at = 1S 

not below the dashed line of G^ WR = 50 dB in Figure 2.10, since the 

resulting ^ does not satisfy Equation (2.25). Thus, the i^G^ value 

at = 0.3 is obtained from Figure 2.9 and is 320, above the dashed 

5 

line of G^ R = 50 dB. Since G^ R = 10 (50 dB), the calculated value 

of ^3 is 3.2°, which also satisfies Equation (2.24). When satellites 
SW and SI are near zenith, and earth station EWR is near the equator, 
the conversion of the topocentric angle fg into the geocentric angle may 
be necessary by use of Equation (2.21). 

H. PLOTS INVOLVING LONGITUDE AND LATITUDE 

The minimum allowable satellite spacings in Figure 2.6 are shown as 
a function of the off-axis angle This angle is the natural 
coordinate for the calculation and allows the construction of universal 
curves, but it may not be optimum for system design. A presentation 
using earth station longitudes and latitudes directly can-also be 
helpful. As an example, the results for a universal system parameter 
value Ron = 35 dB are replotted in Figures 2.11, 2.12, 2.13, and 2.14, 
in which the longitude of SW relative to EWR is used as one coordinate 
and the required satellite separation is the other. Note that the 
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satellite separations >are' shown in terms of the geocentric angle. 
However, such a 'presentation* is no longer universal. As shown in Figure 
2.2, EWR is located inside the coverage area of SW and is not 
necessarily at the ground aim point of SW, but EIR designates the ground 
aim point of SI, whether or not an earth station is actually located 
there. The C/I ratio at EWR is used to determine the allowable 
satellite separations. In Figures 2.11 and 2.13, the longitude of EIR 
relative to EWR is shown as an explicit parameter, while the antenna 
gains, the earth station latitudes, and the universal system parameter 
Rqn are implicit parameters in the calculation. In Figures 2.12 and 
2.14, the latitude of EIR is shown as an explicit parameter, while the 
antenna gains, the latitude of EWR, the longitude of EIR and the 
universal system parameter Rpfj are implicit parameters in the 
calculation. The parameter values Rjjfj = 35 dB, G$n = 50 dB, and G^^R = 
50 dB are used in Figures 2.11 and 2.12, and Rqm = 35 dB, Gsij = 40 dB, 
and G^r = 50 dB are used in Figures 2.13 and 2.14. 

The calculation method regarding these four figures is based on 
Equation (2.14). First, the locations of EWR and EIR and the parameter 
values of Rq^, G$it, and Gewr are fixed. Then, for each location of SW 
which is visible from EWR, the locations of SI (both east and west 
relative to SW) which satisfy Equation (2.14) are determined by means of 
an iterative search. Also, SI must be visible from both EWR and EIR in 
this calculation. The computer program used to produce Figures 2.11, 
2.12, 2.13, and 2.14 is given in Appendix E. 
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Consider first the case of the earth stations being collocated 
(satellite antenna beams being overlapped). For the assumed system 
parameters, the required geocentric satellite separation angles for EWR 
at 0°N, 20°N, 40°N, and 60°N are 4.1°, 4.2°, 4.4°, and 4.6°, 
respectively, at the highest elevation of SW as seen from its earth 
station (relative longitude of SW 0°E), and they all increase slightly 
to 4.8° at the lowest elevation of SW, i.e., these curves are fairly 
flat, and the elevation of SW has very little effect on the required 
satellite spacing when the antenna beam of SI is not restricted. 
Therefore, it can be said that, when satellites which use all allowed 
frequency bands and have no service area (beam) restrictions are placed 
over the orbital arcs with constant separation (about 4.8° in this 
example), the orbit would be "full". Also, note that since ip2 1S a 1 wa Y s 
equal to 0° in this case, the corresponding topocentric satellite 
spacing angle obtained from Figure 2.6 or Figure 2.9 is 4.8°, regardless 
of the satellite elevation. This illustrates that the difference 
between the geocentric angle and the topocentric angle increases 
slightly with increasing satellite elevation. 

Consider next the case of the earth stations being separated 
(satellite antenna beams being restricted). It is evident that much 
closer satellite spacing along the orbital arc is achievable with 
increasing the service area separation. Figures 2.11 and 2.13 show that 
the service area separation in the longitudinal direction greatly 
reduces the required satellite spacing at moderately high elevations of 
SW, whereas Figures 2.12 and 2.14 show that the service area separation 
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Figure 2.11 Minimum geocentric satellite spacing angles when earth stations are separated 
in longitudinal direction. = 35 dB, Gsjj = GfryR = 50 dB, with EWR and EIR 
at (a) 0°N of latitude, (b) 20°N, (c) 40°N, and (d) 60°N. 
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Figure 2.11 (Continued) 
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' LONGITUDE OF SW RELATIVE TO EUR LONGITUDE OF SW 


LATITUDE OF EIR j “ LATITUDE OF EIR 



Figure 2.14 Minimum geocentric satellite spacing angles when earth stations are separated 
in latitudinal direction. Rdn = 35 dB, G$it = 40 dB, and G^r = 50 dB, with 
EIR at the same longitude as EWR and EVJR at (a) 0°N of latitude (b) 20°N, (c) 
40°, and (d) 60°N. 





Figure 2.14f (Continued) 



in the latitudinal direction also reduces the required satellite spacing 
but is not greatly affected by the elevation of SW except in the case of 
service areas located at high latitudes. Note that the service area 
separation implies the size of the off-axis angle ^2 rather than the 
actual distance between the earth stations, and that the required 
satellite spacing decreases with increasing ^ 2 * Therefore, it is 
important to consider not only an adequate geographical separation 
between the service areas but also the satellite elevations in order to 
make the off-axis angle ^2 as large as possible and hence reduce the 
required satellite spacing. 

J. EARTH STATION ANTENNA-PATTERN EFFECTS ON SATELLITE SPACING 

When the main beams of SI and SW overlap each other, or when SI is 
intended to serve a very wide area such as earth coverage, the 
allowable satellite spacing is determined only by the receiving antenna 
of EWR, i.e.. Equation (2.14) is reduced to 

R DN = ‘ D EWR^3’ G EWR^ dB • ( 2 » 26 ) 

For the Ron values of interest, only the sidelobe region of the earth 
station antenna designated by (E-3) makes a contribution in Equation 
(2.26), i.e., 

R DN = G EWR “ 32 + 25 1 og 1 0 ( ^ 3 ) dB , (2.27) 

By assigning another universal system parameter as 
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(2.28) 


Y DN " " R DN + G EWR dB 


- “(C/ 1 ) EWR + e swT " E SIT + G EWR dB 


Equation (2.27) is rewritten as 


Y DN “ 32 ~ 25 logio(ip-j) dB 


(2.29) 


Suppose earth station antennas with a better sidelobe level are 
available, so that 

G EWR^ ^3* G EWR^ = 29 - 25 logm^j) ” G EWR dB (2.30) 

which is given in CCIR Rec. 580 [10], Then Equation (2.29) is replaced 
by 

Y dn = 29 - 25 log 10 ( <|> 3 ) dB . (2.31) 


Figure 2.15 shows the graphical results based on Equations (2.29) and 
(2.31). It is presented in the same format as in Figures 2.11 through 
2,14 except that Yq^ is used as a parameter instead of Ron, and there is 
therefore no need to specify the antenna gain values. The computer code 
in Appendix E is used to produce Figure 2.15. The program, as written, 
produces the solid curves corresponding to Equation (2.29); in order to 
produce the dashed curves corresponding to Equation (2.31), the value 32 

in the subroutine REAL FUNCTION DISC ^E in Appendix E is replaced with 

the value 29. It is obvious that better discrimination in the sidelobe 
region can provide closer satellite spacing. 
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Figure 2.15 Minimum geocentric satellite spacing angles when earth stations are collocated 
(satellite antenna beams are overlapped). Solid lines represent satellite 
spacing angles corresponding to sidelobe gain 32-251 og i o ( ^ 3 ) » while dashed 
lines correspond to 29-251ogio( iM . Y nN i s defined in Equation (2.28). 
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CHAPTER III 


UP-LINK INTERFERENCE 


The single-entry C/I ratio on the up links will now be derived in a 
manner similar to that used for the down links. The interference 
geometry between up -link networks is shown in Figures 3.1 and 3.2 with 
the geometric parameters of interest. The wanted signal is received 
from earth station EWT, while the interference signal arrives from earth 
station EIT. The calculation of the C/I ratio is performed at satellite 
SW. Note that, since there are a multiple number of earth stations 
transmitting to a satellite on the up link, it is not possible for its 
satellite antenna to be pointed directly at all the earth station 
antennas, i.e., the receiving antenna of SW is not necessarily pointed 
directly at the transmitting antenna of EWT. However, the EWT and EIT 
antennas are assumed to be pointed directly at their respective 


satellites as required by good system design. 


The wanted carrier power C$vi and the interference signal power Isw 
received at satellite SW are obtained by application of the Friis 
transmission formula [3] as 



and 


P G 
^EWT^EWT 


{G SWR D SWR^6 ,G SWR^ }A W 
(4irb) ^ 


(3.1) 


l SW 


P EIT {G EIT D EIT^4’ G EIT^ HG SWR D SWR^1’ G SWR^ } x 
( 4ird) ^ 


I 


(3.2) 
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Figure 3.1 Interference geometry between up-link networks. 
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where the symbols, notation, and units adhere to the same convention as 
in the down -link case. 

The satellite antenna reference pattern of Figure 2.3 is used to 
determine the value of D$wr( 'PIjGswr) with appropriate parameter 
changes: ^ and Gswr on the up link correspond to ^2 anc ' ^SIT on t * ie 

down link, respectively. Similarly, the earth station antenna 
reference pattern of Figure 2.4 is used to determine the value of 
d EIT( M»^EIT)» with ^4 and on the up link corresponding to ^3 and 
Gewr on the down link, respectively. 

Dividing Equation (3.1) by Equation (3.2) gives the C/I ratio at 
satellite SW as 


(C/I) 


p ewt g ewt 


°swr ( %» g swr ) 


SW 


P EIT G EIT ^WR^I’SjR^EIT^ V G EIt) b 2 


(3.3) 


p p 

Although the factor c /a was eliminated from the down-link equation by 

assuming c^/a 2 = 1 , the correspondi ng factor d 2 /b 2 cannot be dropped 

from Equation (3.3) because even for small values of the difference 

between the distances b and d is not necessarily negligible, and the 
2 2 

assumption of d /b = 1 is not appropriate. Instead, since the quantity 
2 2 

d /f is in the range of 1 t 0.028 for satellite spacing angles of 5° or 

2 p 

less, the factor d can be replaced by f , and then Equation (3.3) is 


rewritten as 


(C/I) 


SW- = 


P SW D SWr(V G SWr) 

P SI °SWR ( V G SWR )D EIT ( V^If) 


(3.4) 
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where 


P SW = 


P G 
^EWTEWT 

4nb^ 


(3.5) 


is the power density at SW due to EWT, and 
P EIT G EIT 

p SI =-^- (3.6) 

is the power density at SI due to EIT. Defining a universal system 
parameter on the up links as 

P SI 1 

R|ip = (C/Dcu n D ( G ) (3.7) 

uk bw p sw u SWR iv b swR ; 

simplifies Equation (3.4) to 


1 

R U p D SWR^1» G SWR )D EIT (,, '4» G EIT^ 


(3.8) 


As in the down-link case, DsWR( ^ j^SWR) in Equation (3.7) can be 
estimated between 0 and -3 dB, depending on how closely EWT is located 
to the SW ground aim point. 

A few words may be in order with regard to the universal system 
parameter Ryp. On the up link, both the wanted and the interfering 
signals, which are transmitted by earth stations located at different 
points on the Earth, will fluctuate due to local precipitation but will 
not, in general, be attenuated to the same degree. Thus, attenuation 
factors which take the local precipitation effect into account should be 
introduced in calculating the C/I ratio for the up link. Such factors 
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are to be included in the power densities P$w and Psi (although not 
included in Equations (3.5) and (3.6)). The universal system parameter 
Ryp, which is represented in terms of the power densities P$w and Psi* 
also includes the local precipitation effect. The relationship between 
rain rate and attenuation is discussed in Appendix F. On the down link, 
both the wanted and the interfering signals travel through approximately 
the same disturbed region. Therefore, their respective attenuations are 
likely to be correlated, and the local precipitation effect may usually 
be neglected. 

Compare the up link-equation and the down-link equation shown 
below, 

R UP = " D SWR^1» G SWR^ " °EIT^4 ,G EIT^ dB ( 3 * 9 ^ 

and 

R DN = ‘ D SIT^2 ,G EIT^ " D EWR^3» G EWR^ dB * (3.10) 

It is evident that these two equations are of exactly the same form. 

This implies that the up-link and the down-link problems can be treated 
as dual problems, and that the contour curves produced for the down link 
can also be used for the up link by simply changing the down-link 
parameters R QN , lj^, i|> 3 , G SIT , and G £WR to the up-link parameters R up , 

^1* ^4* G SWR* and g e I T * respectively. Thus, all the plots generated for, 
the down link in Chapter II can also apply to the up link with the 
appropriate parameter changes. 

The universal curves for the minimum allowable satellite spacing 
based on up-link considerations are shown in Figures 3.3 and 3.4. The 
values of obtained from Figures 3.3 and 3.4 must satisfy 
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(3.11) 


-U.J- 

i|> > 26.3 6 degrees 

A EIT . 

and 

< 20 y 5.35 + 5 log io(d/ X) / (d/X) degrees, (3.12) 

respectively. For convenience, the dashed lines in the figures indicate 
where to switch from one to the other in the particular case of G^jj = 

50 dB and 60 dB. Figure 3.3 is valid above the dashed lines for the 
given values of G^u, and Figure 3.4 is valid below the dashed lines. 

The following observation can be made regarding the up link. When 
the transmitting earth station EIT is located at the ground aim point 
of the SW receiving antenna, i.e., = 0°, the worst satellite spacing 

results. Separating EIT from the ground aim point of SW, that is, 
increasing the off-axis angle fj, reduces the required satellite 
spacing. The restriction on the transmitting earth station locations, 
i.e., the service area assignments on the up link, turns out to be 
significant in allowing closer satellite spacing. 
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Figure 3.3 First part of universal curves for minimum allowable 
satellite spacing angle ty. as function of normalized 
off-axis angle ^/^g, based on up-link considerations. 

^4 obtained from this figure must satisfy Equation (3.11) # 
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Figure 3.4 Second part of universal curves for minimum allowable 
satellite spacing angle as function of normalized 
off-axis angle based on up-link considerations. 

i |>4 obtained frornthis figure must satisfy Equation (3.12). 
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CHAPTER IV 


CONCLUSION 


This study deals with the minimum allowable satellite spacing for 
FSS systems which attains an acceptable single-entry C/I ratio, 
using co-pol ari zed circular antenna beams conforming to Recommendations 
and Reports of the CCIR. 

Universal curves for the minimum allowable satellite spacings are 
constructed, in which the natural coordinates are normalized off-axis 
angles, defined as the angle between a main beam axis and an 
interference direction, seen at the antenna. In the down-link case, the 
independent (input) off-axis angle is that of the closest wanted earth 
station, seen at the interfering satellite; in the up-link case, it is 
that of the closest interfering earth station, seen at the wanted 
satellite. The dependent (output) angle in each case is the topocentric 
satellite spacing angle; however, the difference between the topocentric 
angle and the geocentric angle is small, the topocentric angle is always 
greater than the geocentric angle, and hence the topocentric angle used 
as the minimum satellite spacing is conservative in this sense. The key 
to the universal representation is the combination of the 
multiplicative/divisor system parameters into a universal factor, which 
can be set equal to the total antenna discrimination required. The 
earth station separation then determines the part of this discrimination 
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supplied by a satellite antenna; the remainder must be supplied by an 
earth station antenna, and this leads to the required satellite 
separation. 

Plots involving earth station longitudes and latitudes directly are 
also presented for some particular parameter values. When the satellite 
beams serve very wide areas such as earth coverage, the allowable 
satellite separation is achieved mainly by the earth station antenna 
discrimination patterns and is almost independent of the elevation 
angle. However, the allowable satellite spacing can be decreased 
substantially by defining service areas for the satellites; thus the 
assignment of limited service area turns out to be an important concept 
for the effective use of the geostationary satellite orbit. For 
longitudinal service area separations, the allowable satellite spacing 
is decreased most when the satellites are near zenith; for latitudinal 
service area separations, the improvement is relatively independent of 
the satellite location except in the case of service areas at high 
latitudes, which is of little practical importance. 

For the case of no service area separation, universal curves are 
constructed to show the satellite spacing improvement available through 
the improvement of the earth-station sidelobe level. 

In general, the considerations of both down-link and up-link 
interference are necessary to determine orbital satellite locations for 
the FSS systems. It is shown that the up-link and down-link 
calculations are dual, so that the same curves and computer codes can be 
used when suitable parameter substitutions are made. 
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APPENDIX A 


PROGRAM KAZU.FOR 

This program computes minimum topocentric satellite spacing 
angle (PSI3) as a function of off-axis angle (PSI2). 


C 

C 

10 

20 

$ 


C 


c 


200 

C 

30 


40 


300 

100 


REAL GAINSI, GAINEW 

REAL RJSIVEN, R_TEST 

REAL PSI 2, PSI 3_C0ARSE , PSI3_FINE 

INTEGER I, J, K 


WRITE(6,10) 

FORMAT (// 1 Enter R in dB.‘) \> 

READ(5,*) R GIVEN . 

WRITE(6,20) 

F0RMAT(//‘ Enter satellite (SI) antenna gain and 1 
/' earth station (EWR) antenna gain in dB . ' ) 
READ( 5,*) GAINSI, GAINEW 



0PEN(UNIT = 7, FILE = 'KAZU.DAT', TYPE = 'NEW') 
DO 100 I = 1, 200 
PSI2 = 0.05 * (I - 1) 


DO 200 J = 1, 100 
PSI3_C0ARSE = 0.1 * J 

R TEST = -DISC S(PSI2, GAINSI) - DISC_E(PSI3_C0ARSE, GAINEW) 
IF(R_TEST .GT.“fc_GIVEN) GO TO 30 
CONTINUE 
GO TO 100 


DO 300 K = 1, 101 

PSI3 FINE = PSI3 COARSE - 0.1 + (K - 1) * 0.001 
R TEST = -DISC STPSI2, GAINSI) - DISC E(PSI3_FINE, GAINEW) 
IF(R TEST .GE. R GIVEN) THEN 
WRTTE(6,40) PSI 2, PS 1 3 FINE 
WRITE(7,40) PSI2, PS 1 I NE 
F0RMAT(2F10.2) 

GO TO 100 
END IF 
CONTINUE 
CONTINUE 
CLOSE ( UNIT = 7) 

STOP 

END 
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Function DISC_S calculates relative antenna gain of a satellite 
in dB from off-axis angle ( PS I ) and antenna gain (GAIN). 

REAL FUNCTION DISC S(PSI.GAIN) 

HPBW = SQRT( 27000. / 10.** (GAIN/ 10. ) ) 

PSI HP = PSI / HPBW 
IF(FSI HP .LE. 1.291) 

$ DISC S = -12. * PSI HP**2 

IF(PSI"HP .GT. 1.291 .AND. PSI_HP .LE. 3.1623) 

$ DISC S = -20. 

I F ( PS I HP .GT. 3.1623) 

$ DISC - S = -7.5 - 25. * L0G10(PSI_HP) 

IF(DISZT S .LE. -GAIN - 10.) 

$ DISCJS = -GAIN - 10. 

RETURN 

END 

Function DISC_E calculates relative antenna gain of an earth 
station in dB from off-axis angle (PSI) and antenna gain (GAIN). 

REAL FUNCTION DISC E(PSI,GAIN) 

D W = SQRT(10.**(G7TIN/10.) / 5.428) 

P3T M = 20. * SQRT(GAIN - 2. - 15. * L0G10(D_W)) / (D_W) 

PSI - R = 15.85 * D W** (—0.6) 

IF(FSI .LT. PSIJJ 

$ DISC E = -2.5 E-3 * (D W * PSI )**2 
IF(PSI“.GE. PSI M .AND. FSI .LT. PSI R) 

$ DISC E = 2. +15. * L0G10(D W) - GFlN 
IF (PSI .GE.PSI R .AND. PSI .LT. 48.) 

$ DISC_E = 32.“- 25. * LOGIO(PSI) - GAIN 
IF(PSI .GE. 48.) 

$ DISC E = -GAIN - iO. 

RETURN 

END 
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APPENDIX B 


A conversion formula from the topocentric angle ^ to the 

geocentric angle <j> g is derived as follows. Referring to Figure B.l, 

/ 

sinU g /2) = U/2)/R . (B.l) 

For small <j>g, 

( 8 . 2 ) 

Since a » c for small 

a c (B. 3) 

*g «^Ty2 ( 1-cos a ^V 2 ( 1-cos * t ) . 




FIGURE B.l Relationship between topocentric angle and geocentric 
angle. 
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for smal 1 t , 


Since cos^ 


1 - 


% 


Y 2 ( 1-cos ip t ) a 


(B.4) 


Therefore, 


% * R \ * R \ » 


(B.5) 


and 


if -1 V 1+ (0 2 -r cos( ^ )co ^ 


= ^ 1.023 - 0.302 cosffy,-^) cose E , 


(B.6) 


where ^ is a midpoint (longitude) between two satellites, and 0£ and <j> E 
are latitude and longitude of an observation point on the Earth, 
respectively. 
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APPENDIX C 


C PROGRAM UNIV1.FOR 

C This program computes PSI3*(GAIN**0.4) as a function of PSI2/HPBW 
C 

REAL R_GIVEN, R_TEST 
REAL PS 1 2 HP, ARG COARSE, ARG_FINE 
C 

INTEGER I, J, K 
C 

WRITE ( 6 , 10) 

10 FORMAT (// 1 Enter R in dB.') 

READ( 5,*) R GIVEN 
C 

OPEN(UNIT = 7, FILE = 'UNIV1.DAT', TYPE = 'NEW') 

C 

DO 100 I = 1, 101 
PSI2 HP = 0.05 * (I - 1) 

C 

DO 200 J = 1, 800 
ARG COARSE = 1. * J 

R TEST = -DISC E(ARG COARSE) - DISC $(PSI2 HP) 

TF(R TEST .GT.“R GIVEN) GO TO 30 
200 CONTINUE “ 

GO TO 100 
C 

30 DO 300 K = 1, 101 

ARG FINE = ARG COARSE - 1. + (K - 1) * 0.01 
R TEST = -DISC - E{ARG FINE) - DISC S(PSI2 HP) 

IF(R_TEST .GE.“R_GI\0) THEN 
WRITE{6,40) PS 1 2 HP, ARG FINE 
WRITE(7,40) PSI2HP, ARG"EINE 
40 FORMAT (2F10.2) “ “ 

GO TO 100 
END IF 
300 CONTINUE 

100 CONTINUE 

CLOSE (UNIT = 7) 

STOP 

END 
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C Function DISC S computes relative gain of a satellite antenna 
C in dB _f rom normalized off-axis angle (PSI/HPBW). 

C 

REAL FUNCTION DISC S(PSI_HP) 

IF ( PS I HP .LE. 1.251) 

$ Disrs = -12. * PSI HP**2 

I F ( PS I TTP .GT. 1.291 MU. PSI_HP .LE. 3.1623) 

$ DISC S = -20. 

IF (PSI TTP .GT. 3.1623) 

$ DISC S = -7.5 - 25. * LOG 10 ( PSI HP ) 

RETURN 
END 

Function DISC E computes relative gain of second si delobe of 
an earth station antenna. 

REAL FUNCTION DISC E(ARG) 

C ARG = PS 1 3 * (GAIN T *0.4) 

DISC E = 32. - 25. * LOGIO(ARG) 

RETURN 

END 



APPENDIX D 


C PROGRAM UNI V2. FOR 

C This program computes PSI3*{d/1ambda) as a function of PS12/HPBW 
C 

REAL R GIVEN, R TEST 
REAL PSI2_HP, AiTG_COARSE, ARG_FINE 
C 

INTEGER I, J, K 
C 

WRITE (6, 10) 

10 FORMAT (// ' Enter R in dB.') 

READ(5,*) R_GIVEN 
C 

OPEN(UNIT = 7, FILE = 'UNIV2.DAT*, TYPE = 'NEW') 

C 

DO 100 1=1, 101 

PSI2 HP = 0.2 * (I - 1) 

C “ 

DO 200 J = 1, 150 
ARG COARSE =1. * J 

R TEST = -DISC_E(ARG COARSE) - DISC S(PSI2 HP) 

IF(R_TEST .GT. RjGIVEN) GO TO 30 " 

200 CONTINUE 

GO TO 100 
C 

30 DO 300 K = 1, 101 

ARG FINE = ARG_COARSE - 1. + (K - 1) * 0.01 
R TEST = -DISC E(ARG FINE) - DISC S(PSI2 HP) 

IE(R_TEST .GE. R_GIVEN) THEN 
WRITE( 6,40) PSI2_HP, ARG FINE 
WRITE(7,40) PS 1 2 HP, ARG~FINE 
40 FORMAT (2F 10. 2) ” 

GO TO 100 - 

END IF 
300 CONTINUE 

100 CONTINUE 

CLOSE (UNIT =7) 

“ " STOP 

END 
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Function DISC S computes relative gain of a satellite antenna 
in dB from normalized off-axis angle (PSI/HPBW). 

REAL FUNCTION DISC_S(PSI_HP) 

IF(PSI HP .LE. 1.291) 

$ DISC S = -12. * PSI HP**2 

I F ( PS I JTP .GT. 1.291 JND. PSI_HP .LE. 3.1623) 

$ DISC S = -20. 

I F ( PS I HP .GT. 3.1623) 

$ DISC S = -7.5 - 25. * L0G10(PSI_HP) 

RETURN 
END 

Function DISC_E computes relative gain of main lobe of an earth 
station antenna. 

REAL FUNCTION DISC E(ARG) 

C ARG = PS 1 3 * (d/lambda) 

DISC_E = -2.5 E-3 * (ARG)**2 

RETURN 

END 



APPENDIX E 


C PROGRAM MURA. FOR 

C This program computes minimum geocentric satellite angles 
C when locations of two earth stations are fixed. 

C 

REAL RGSO RE 

REAL LONGSW, LONGS I COARSE, LONGSI FINE 
REAL LONGEW, LATIEW, LONGEI, LATIET 
REAL GAINS I, GAINEW 
REAL R GIVEN, R TEST 

REAL L EWEI , L TWEW, L_SIEW, LJSIEI, L_SWSI 
REAL PTI2, PS IT 
REAL EL EWSI , EL EISI 
REAL RATTGE1, RANTE2, SEPA, DELTA 
C 

INTEGER I, J, K, L, M, N 
C 

RGSO = 4.2152 E7 
RE = 6.371 E6 
C 

WRITE(6,10) 

10 FORMAT(//' Enter earth station (EWR) location in degrees.', 

$ /' (longitude, latitude)') 

WRITE(6,20) 

20 F0RMAT('0', ' The longitude value is positive for the eastern', 

$ 'hemisphere', /' and negative for the western', 

$ ' hemisphere. Likewise,', /' the latitude value is', 

$ ' positive for the northern hemisphere' /' and negative', 

$ ' for the southern hemisphere.') 

READ(5,*) LONGEW, LATIEW 
WRITE(6,30) 

30 F0RMAT(//' Enter earth station (EIR) location in degrees.', 

$ /' (longitude, latitude)') ------ - • 

READ( 5,*) LONGEI, LATIEI 
WRITE(6,40) 

40 F0RMAT(// ' Enter R in dB.') 

READ(5,*) R GIVEN 
WRITE(6, 50) 

5Q F0RMAT(//' Enter satellite (SI) antenna gain and', 

$ /' earth station (EWR) antenna gain in dB.') 

READ(5,*) GAINSI , GAINEW 
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c 

0PEN(UNIT = 7, FILE = 'MURAE.DAT', TYPE = 'NEW') 

0PEN(UNIT = 8, FILE = 'MURAW.DAT', TYPE = 'NEW') 

Satellite (SW) can be seen from earth station (EWR) between 
(RANGE 1 ) deg. and (RANGE2) deg. of longitude. 

RANGE 1 = LONGEW - ACOSD(RE/(RGSO*COSD(LATIEW) ) ) 

RANGE 2 = LONGEW + AC0SD(RE/ (RGSO*COSD(LATIEW) ) ) 

DELTA = (RANGE2 - RANGE 1) / 200. 

L_EWE I = DIST(RE, LATIEW, LONGEW, RE, LATIEI, LONGEI ) 

The case of satellite (SI) being located in the east of 
satellite (SW). 

DO 100 I = 1, 201 
L'ONGSW = RANGE 1 + DELTA * (I - 1) 

C 

L SWEW = DIST (RGSO, 0., LONGSW, RE, LATIEW, LONGEW) 

C 

DO 200 J = 1, 20 
LONGS I COARSE = LONGSW + J * 0. 5 
C 

L_SIEW = DIST(RGS0, 0., LONGS I_C0ARSE , RE, LATIEW, LONGEW) 
L_S I E I = DIST (RGSO, 0., LONGSI COARSE, RE, LATIEI, LONGEI) 
L SWSI = DIST(RGS0, 0., L0NGSW7 RGSO, 0., LONGSI COARSE) 

C 

PS 1 2 = OFF AX I S ( L SIEW, L SIEI, L EWEI) 

PSI3 = OFF AXIS (L — SWEW, L"SIEW, L SWSI) 

C _ _ _ _ 

R_TEST = -DISC_S(PSI2,GAINSI ) - DISC_E (PSI3,GAINEW) 

I F ( R TEST .GT. R GIVEN) GO TO 60 
200 CONTINUE 
GO TO 100 
C 

60 DO 300 K = 1, 101 

LONGSI FINE = LONGSI COARSE - 0.5 + (K - 1) * 0.005 
C “ 

L_S I EW = DIST(RGS0, 0., LONGSI FINE, RE, LATIEW, LONGEW) 
LJSIEI = DIST (RGSO, 0., LONGSI_FINE, RE, LATIEI, LONGEI) 

L SWSI = DIST(RGSO, 0., LONGSW, RGSO, 0., LONGSI_FINE) 

C 

PS 1 2 = OFF AXIS(L SIEW, L SIEI , L_EWEI ) 

PS 1 3 = OFF AXIS(L SWEW, L SIEW, L SWSI) 

C “ 

R TEST = -DISC S(PSI 2.GAINSI ) - DISC E(PSI3,GAINEW) 

C ” 

I F ( R TEST .GE. R GIVEN) GO TO 70 
300 CONTINUE 
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70 EL EWSI = ELEV(LATIEW, LONGEW, LONGSI FINE) 

EL - EISI = ELEV(LATIEI , LONGEI, LONGSQINE) 

C 

IF(EL EWSI .GE. 0. .AND. EL EISI .GE. 0.) THEN 
SEP7T = LONGSI FINE - LONGSW 
WRITE(6,80) LlJNGSW, SEPA 
WRITE(7,80) LONGSW, SEPA 

80 F0RMAT(2F10.2) 

END IF 

100 CONTINUE 

The case of satellite (SI) being located in the west of 
satellite (SW). 

DO 400 L = 1, 201 
LONGSW = RANGE 1 + DELTA * (L - 1) 

C 

L_SWEW = DIST(RGSO, 0., LONGSW, RE, LATIEW, LONGEW) 

C 

DO 500 M = 1, 20 
LONGSI COARSE = LONGSW - M * 0.5 
C 

L SIEW = DIST(RGSO, 0., LONGSI COARSE, RE, LATIEW, LONGEW 
LJ5IEI = DIST(RGSO, 0., LONGS I~COARSE, RE, LATIEI, LONGEI 
L SWSI = DIST(RGSO, 0., LONGSW, RGSO, 0., LONGS I_COARSE) 

C 

PS 1 2 = OFF AXIS(L SIEW, L SIEI, L EWEI) 

PS 1 3 = OF F~A X I S ( L~ SWE W , L SIEW, L"SWSI) 

C “ 

R TEST = -DISC S(PSI2,GAINSI ) - DISC E(PSI3,GAINEW) 

C ” 

IF (R TEST .GT. R GIVEN) GO TO 90 
500 CONTINUE 

GO TO 400 
C 

90 DO 600 N = 1, 101 

LONGSI FINE = LONGSI COARSE + 0.5 - (N - 1) * 0.005 
C 

LjSIEW = DIST(RGSO, 0., LONGSI FINE, RE, LATIEW, LONGEW) 

L SIEI = DIST (RGSO, 0., LONGSI~FlNE, RE, LATIEI, LONGEI) 
L“SWSI = DIST(RGSO, 0., L0NGSW7 RGSO, 0., LONGSI_FINE) 

C 

PSI2 = OFF AXIS(L SIEW, L SIEI, L EWEI) 

PSI3 = OFF - AX I S ( L — SWE W , L SLEW, L“SWSI) 

C " 

R TEST = -DISC S(PSI2,GAINSI) - DISC E(PSI3,GAINEW) 

C 

IF(R_TEST .GE. RjGIVEN) GO TO 95 
600 CONTINUE 
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95 EL EWSI = ELEV(LATIEW,. LONGEW, LONGS I FINE) 

EL“EISI = ELE V ( LAT I E I , LONGEI , LONGSIlFINE) 

C 

I F ( EL EWSI .GE. 0. .AND. EL EISI .GE. 0.) THEN 
SEP7T = LONGS I FINE - L0NG5W 
WRITE(6,80) LUNGSW , SEPA 
WRITE(8,80) LONGSW, SEPA 
END IF 

400 CONTINUE 

CLOSE (UNIT = 7) 

CLOSE (UNIT = 8) 

STOP 
END 

Function DISC_S computes relative gain of a satellite antenna 
in dB from off-axis angle ( PS I ) and antenna gain (GAIN). 

REAL FUNCTION DISC S(PSI,GAIN) 

HPBW = SQRT (27000. / 10.**(GAIN/10.) ) 

PSI HP = PSI / HPBW 
I F ( PS I HP .LE. 1.291) 

nisr s = -i? * psi hp**2 
I F ( PS I HP .GT. 1.291 .BID. PSI HP .LE. 3.1623) 

$ Disr_s = -20. 

IF ( PSI HP .GT. 3.1623) 

$ DISC S = -7.5 - 25. * L0G10(PSI_HP) 

IF(DISC“S .LE. -GAIN - 10.) 

$ DISC S = -GAIN - 10. 

RETURN ” 

END 

"Ftincti orPDlSC r E^comput^i^fel'ati^ ve — gai rroTluF^ 
antenna in dB Trom off-axis angle (PSI) and antenna gain (GAIN) 

REAL FUNCTION DISC E(PSI,GAIN) 

D_W = SQRT( 10.**(G7TIN/10. ) / 5.428) 

PSI M = 20. * SQRT (GAIN - 2. - 15. * L0G10(D_W)) / D_W 
PSI R = 15.85 * D W**(-0.6) 

I F (PS I .LT. PSI MT 

$ DISC E = —2 . 6 E— 3 * (D W * PSI )**2 
I F ( PS I . GE . PSI M .AND. HSI .LT. PSI R) 

$ DISC_E = 2. +”15. * L0G10(D_W) - G3HN 
I F ( PS I .GE. PSI R .AND. PSI .LT. 48.) 

$ DISC_E = 32. “ 25. * LOGIO(PSI) - GAIN 
I F ( PS I .GE. 48.) 

$ DISC E = -GAIN - 10. 

RETURN - 
END 
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Function DIST computes distance between two points 
represented in the cylindrical coordinate system. 

(radius, latitude, longitude) 

REAL FUNCTION DIST(R1, THETA1, PHI1, R2, THETA2, PHI2) 

XI = R1 * COSD(THETAl) * COSD(PHIl) 

Yl = R1 * COSO(THETAl) * SIND(PHIl) 

Z1 = R1 * SIND(THETAl) 

X2 = R2 * C0SD(THETA2) * C0SD(PHI2) 

Y2 = R2 * C0SD(THETA2) * SIND(PHI2) 

Z2 = R2 * SIND(THETA2) 

DIST = SQRT ( ( Xl-X2)**2 + (Yl-Y2)**2 + (Zl-Z2)**2) 

RETURN 

END 

Function ELEV computes elevation angle of a satellite (PHI_S) 
seen from an earth station (THETA_E, PHI E ) . 

REAL FUNCTION ELEV(THETA_E, PHI_E, PHI_S) 

RGSO = 4.2152 E7 
RE = 6.371 E6 

IF(SQRT (1. - (COSD (THETA E))**2 * (C0SD(PH1 E-PHI_S) )**2) 

$ .EQ. 0.) THEN 
ELEV = 90. 

ELSE 

ELEV = ATAND(( COSD (THETA E) * C0SD(PHI E-PHI S) - (RE/RGSO) ) 
$ / (SQRT (1. - (COSTT(THETA E))**2“* (C030(PHI_E-PHI_S) ) 

$ ** 2 ) ) ) 

END IF 
RETURN 
END 

Function 0FF_AXIS computes off-axis angle ( PS I ) . 

REAL FUNCTION OFF AXIS(X, Y, Z) 

ARG OFF AXIS = (X**2 + Y**2 - Z**2) / (2. * X * Y) 

IF(ARG TTFF_AXIS .LT. 1. ) THEN 
OFFJTXIS = AC0SD(ARG OFF AXIS) 

ELSE “ “ . . . - 

OFF AXIS = 0. 

END IF 
RETURN 
END 
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The relationship between rain rate R (mm/hr. )» as measured at the 
surface of the Earth, and specific attenuation a (dB/km) can be 
approximated by the equation [11] 

a = aR b dB/km , (F.l) 

where a and b are frequency- and slightly temperature-dependent 
constants. Table F.l presents a listing of the a and b coefficients for 
several frequencies of interest for satellite communications. Also, 
values of specific attenuation at each frequency are listed for rain 
rates of 10, 50, and 100 mm/hr. For example, assuming an effective path 
length of 4 km through rain, the magnitude of the rain attenuation is 
1.2 dB at 6 GHz for a 50 mm/hr. rain and is 2.63 dB at 6 GHz for a 100 
mm/hr. rain. Therefore, for regions of bad rain climate, it may be 
necessary to assume a certain value (dB) of rain attenuation in the 
wanted signal and none in the interfering signal and to adjust the C/I 
value used in Ryp accordingly when calculating the allowable satellite 
separation for the up links. 
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TABLE F.l 


COEFFICIENTS a AND b OF EQUATION (F.l) FOR CALCULATION 
OF RAIN ATTENUATION AT RAIN TEMPERATURE OF 0°C [11] 


Frequency, GHz 

Coefficient 

a, dB/km, for R specified 
in mm/hr 


a 

b 

10 

50 

100 

2 

0.000345 


0.003 

0.011 

0.021 

4 

.00147 

1.016 

.015 

.078 

.158 

6 

.00371 

1.124 

.049 

.30 

.657 

12 

.0215 

1.136 

.29 

1.83 

4.02 

15 

.0368 

1.118 

.48 

2.92 

6.34 

20 

.0719 

1.097 

.90 

5.25 

11.24 

30 

.186 

1.043 

2.05 

11.0 

22.7 

40 

.362 

.972 

3.39 

16.2 

31.8 

94 

1.402 

.744 

7.78 

25.8 

43.1 













